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Abstract. Knowledge regarding the coronal magnetic field is important for the 
understanding of many phenomena, like flares and coronal mass ejections. Because of 
the low plasma beta in the solar corona the coronal magnetic field is often assumed to 
be force-free and we use photospheric vector magnetograph data to extrapolate the 
magnetic field into the corona with the help of a non-linear force-free optimization 
code. Unfortunately the measurements of the photospheric magnetic field contain 
inconsistencies and noise. In particular the transversal components (say Bx and By) 
of current vector magnetographs have their uncertainties. Furthermore the magnetic 
field in the photosphere is not necessary force-free and often not consistent with the 
assumption of a force-free field above. We develop a preprocessing procedure to 
drive the observed non force-free data towards suitable boundary conditions for a 
force-free extrapolation. As a result we get a data set which is as close as possible 
to the measured data and consistent with the force-free assumption. 
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1. Introduction 

The structure of the solar corona out to a few solar radii is domi- 
nated by the magnetic field. Knowledge regarding the coronal mag- 
netic field is therefore important to understand physical processes like 
flares and coronal mass ejections. Unfortunately direct measurements 
of the magnetic field are difficult for the following reasons: The high 
plasma temperature in the corona broadens the line profile orders of 
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magnitudes above the Zeeman splitting. In addition coronal lines are 
optical thin and consequently the line-of-sight integrated character of 
the measurements complicates their interpretation. 

It has been proposed to use magnetic sensitive coronal line ob- 
servations in order to constrain the coronal magnetic field ( see, e.g. 
House (1977) , Arnaud and Newkirk (1987), Judge (1998)). Lin et al. 
(2004) demonstrated that Zeeman-effect measurements can be per- 
formed. However, due to the line-of-sight effect the interpretation of 
these lines is difficult. Also, recent measurements of chromospheric and 
lower coronal magnetic fields were made for a few individual cases, 
e.g. (Solanki et al., 2003; Lagg et al., 2004). Such multi-wavelength 
spectropolarimetric observations have not been done routinely in the 
past but are believed to become much more commonplace, and reliable, 
in the near future. 

As an alternative to these direct measurements, methods have been 
developed to use the magnetic field observed on the photosphere for 
an extrapolation into the corona. The extrapolation is not free from 
assumptions regarding the coronal plasma. It is helpful, that the low 
and middle corona contains a low (3 plasma which allows to neglect the 
plasma pressure in first order and use a force-free magnetic field model. 
A force-free magnetic field is characterized by the electric currents 
parallel to the magnetic field lines, i.e., 



where the coefficient a is constant along the field lines but may vary 
between different field lines. Since (1) due to its intrinsic non-linearity is 
difficult to solve, simplifications of (1) are sometimes used like potential 
fields (q = 0, i.e., no currents; see e.g. (Schmidt, 1964; Semel, 1967; 
Cuperman et al., 1989)) and linear force-free fields (a = const; see e.g. 
(Chiu and Hilton, 1977; Seehafer, 1978; Seehafer, 1982; Alissandrakis, 
1981; Semel, 1988; Demoulin and Priest, 1992)). 

The non-linear force-free case (Sakurai, 1981; Wu et al., 1990; Roume- 
liotis, 1996; Amari et al., 1997; Yan and Sakurai, 2000; Wheatland 
et al., 2000; Wheatland, 2004; Regnier et al., 2002; Wiegelmann and 
Neukirch, 2003; Wiegelmann, 2004; Valori et al., 2005) is challenging 
both theoretically due to the non-linearity of the underlying math- 
ematical problem and observationally because a measurement of the 
full magnetic vector on the photospheric is required. For investigations 
of instabilities like filament eruptions, it is essential to consider the 
general non-linear force-free field case, because with the simplified field 
models mentioned above the free magnetic energy which might drive 
the instability cannot be described 



V x B = aB 



V B = 



(1) 
(2) 
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Moreover, a comparison between extrapolations and measured fields 
(Solanki et al., 2003) in a newly developed active region revealed that 
a non-linear model agrees better with the observations than a potential 
and linear force- free magnetic field model (Wiegelmann et al., 2005). 

A non-linear force-free extrapolation of coronal magnetic fields from 
measured photospheric data is a challenging problem for several rea- 
sons. The magnetic field in the photosphere is not necessary force- 
free (see Gary (2001)) and one would rather prefer to use the vector- 
magnetic field at the basic of the corona as boundary condition for a 
non-linear force-free extrapolation, but here measurements are usually 
not available. 

Another problem is the 180° ambiguity in the transversal compo- 
nents of the photospheric data, which has to be removed with e.g. the 
minimum energy method (Metcalf, 1994). An additional complication 
is that the transversal components are more difficult to measure and 
bear much more noise than the line-of-sight component. For the data 
of the Solar-Flare Telescope (SFT) used in section 5, for example, the 
noise level is about 10G for the line-of-sight magnetic field and 100G 
for the transverse field. The full disk vector magnetograph SOLIS (Syn- 
optic Optical Long-term Investigations of the Sun, U.S. National Solar 
Observatory, Kitt Peak) is expected to have a noise level of about 
1G in the line-of-sight and about 50G in the transverse components 
(C. Keller, private communication). The reason for the much higher 
noise in the transversal magnetic field (Bt) compared with the line-of- 
sight field (Bl) can be understood by the way how the photospheric 
magnetic field is derived from the four measured Stokes components. 
Roughly Bl = clV/I where V is the circular polarization and cl is a 
constant. Then 

5B L = cl— oc — oc v 7 !, (3) 

where 81 is the noise in /. This photon noise is independent of Bl- 
Even if the Q and U signals are generally weaker than the V signal, 
the major source of errors in the force balance is the noise in B t : B 2 = 
c2 \JQ 2 + U 2 /I where Q and U are linear polarization intensities and 
c2 is another constant. Then 

9 r ar QSQ±USU n V5Q 2 + 5U 2 SI 

2B t oB t = c2 — — ~ c2 ~ c2 — ~ const. (4) 

VQ 2 + U 2 I I I K ' 

As a consequence the transvers magnetic field noise is in particular 
high in weak field regions. In strong field regions vector magnetic field 
measurements based on full Stokes spectro-polarimetry generally yield 
uncertainties of a few degrees in the orientation angles of the field 
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vector in active regions. This situation is also constantly improving 
with improved instrumentation. 

For the practical computation of the coronal magnetic field it is 
also helpful if the boundary data are smoothed to some degree. Short 
wavelength fluctuations in the surface magnetic field die out rapidly 
with height (for a potential field with exp(— kz) for scales of horizontal 
wavenumber A;). Hence, boundary data with scales much finer than the 
height of the numerical box require a very fine grid to be sufficiently 
resolved but they hardly have a effect on the result, except in a very 
small boundary layer above the surface. To keep the numerical effort 
limited, the boundary data are therefore usually smoothed. In this 
context it should also be noted that the force-free equation (1) is scale 
invariant if a is rescaled appropriately. Hence, by smoothing of the 
boundary data we do not lose any physics which we have not already 
cast away by the restriction to a stationary force-free field model. 

A more fundamental requirement on the boundary data is its con- 
sistency with the force-free field approximation. As has been shown by 
Aly (1989), a balance of the total momentum and angular momentum 
exerted onto the numerical box by the magnetic field leads to a set of 
boundary integral constraints for the magnetic field. These constraints 
should also be satisfied on the solar surface for the field at the coronal 
base in the vicinity of a sufficiently isolated magnetic region and in a 
situation where there is no rapid dynamical development. 

In summary, the boundary data for the force-free extrapolation 
should fulfill the following conditions: 

1. The data should coincide with the photospheric observations within 
measurement errors. 

2. The data should be consistent with the assumption of a force- free 
magnetic field above. 

3. For computational reasons (finite differences) the data should be 
sufficiently smooth. 

Within this work we describe a numerical procedure, written in IDL, 
which pre-processes vector magnetograph data so that the above con- 
ditions are satisfied as close as possible. 

We outline the paper as follows. In section 2 we specify how to 
check if a given measured vector magnetograms is consistent with the 
assumption of a force-free magnetic field above. We describe, how to 
derive consistent force-free boundary conditions for a force-free mag- 
netic field extrapolation from the measured data in section 3. In section 
4 we use a known semi-analytic force-free model to check our method 
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and we apply the method in section 5 to an example of the observed 
vector magnetogram. 

2. Consistency check of vector magnet ograms 

In the following we assume that the 180° ambiguity in the measured 
transverse components has been removed and the magnetogram has 
been observed close to the disk center. "Close" here means that the ver- 
tical component B z in our numerical box can be more or less identified 
with the line-of-sight component and bears a smaller measurement error 
than the horizontal components B x and B y . Since, as stated above, 
the measurement error between the line-of-sight and the transverse 
components differs by about one order of magnitude, the angle between 
line-of-sight and vertical should not exceed ~ 0.1 radian. 

Another a-priori assumption about the photospheric data is that 
the magnetic flux from the photosphere is sufficiently distant from the 
boundaries of the observational domain and the net flux is in balance, 
i.e., 



Generally, the flux balance criterion has to be applied to the whole, 
closed surface of the numerical box. However, we can only measure 
the magnetic field vector on the bottom photospheric boundary and 
the contributions of the lateral and top boundary remain unspecified. 
However, if a major part of the known flux from the bottom bound- 
ary is uncompensated, the final force-free magnetic field solution will 
markedly depend on how the uncompensated flux is distributed over 
the other five boundaries. This would result in a major uncertainty of 
the final force free magnetic field configuration. We therefore demand 
the flux balance to be satisfied with the bottom data alone. If this is 
not the case, we classify the reconstruction problem as not uniquely 
solvable within the given box. 

Aly (1989) used the virial theorem to define which conditions a 
vector magnetogram has to fulfill to be consistent with the assumption 
of a force-free field above. We repeat here the force-free and torque-free 
condition and refer to the paper of (Aly, 1989) for details. 

1. The total force on the boundary vanishes 






B y B z dxdy = 




H2092L_rev.tex; 5/02/2008; 15:07; p. 5 



6 



Wiegelmann et al. 



2. The total torque on the boundary vanishes 

/ x (B 2 + B 2 ) dxdy = / x B 2 dx dy 
JS Js 

f y{B 2 x + B 2 y ) dxdy= f y B 2 Z dx dy 

/ y B X B Z dxdy = \ x B V B Z dx dy 
Js Js 

Note that if condition 1) is fulfilled, the constraints 2) are independent 
on where the origin for x, y is located. In our code the origin is usually 
at the lower left corner of the bottom boundary face. 

As with the flux balance, the Aly integral criteria in general have to 
be applied to the whole surface of the numerical box. Since we assumed 
that the photospheric flux is sufficiently concentrated in the center and 
the net flux is in balance, we can expect the magnetic field on the lateral 
and top boundary to remain small and hence these surfaces will not 
yield a large contribution to the integrals of the Aly-constraints above. 
We therefore impose the Aly-criteria on the bottom boundary alone. 

Aly (1989) already pointed out that the magnetic field is probably 
not force-free in the measured region because the plasma (5 in the 
photosphere is of the order of one and pressure and gravity forces are 
not negligible. We however expect that the observed photospheric field 
is, after removing scales below a super-granular diameter by smoothing, 
representative for the field at the coronal base. In the corona, however 
we have (3 « 1CT 4 and consequently the magnetic field should be close 
to force-free in a stationary situation. 

To quantify the quality of vector magnetograms with respect to the 
above criteria we introduce three dimensionless parameters: 

1. The flux balance parameter 

J s B z dx dy 
£flux = J S \B Z \ dxdy 

2. The force balance parameter ef orce = 

| Js B X B Z dx dy\ + | Js B y B z dx dy\ + | J S (B 2 + B 2 ) - B 2 dx dy\ 
J s {Bl + B 2 + B 2 )dxdy 

3. The torque balance parameter e^ or q Ue = 

| J A x((B 2 x + Bp - B 2 )dxdy\ + | J A y((B 2 x + B 2 ) - B 2 ) dxdy\ + | J s yB x B z - xB y B z dxdy\ 

J s V^+F (B 2 + Bl + Bl) dx dy 
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An observed vector magnetogram is then flux-balanced and consistent 
with the force-free assumption if: £fl ux <C 1, £f orce <C 1, etorque ^ 1- 



3. Method 

Even if we choose a sufficiently flux balanced isolated active region 
( e flux ^ 1) we ^ n< ^ that usuan Y the force-free conditions £f orce <C 
1, ^torque ^ 1 are n °t fulfilled for measured vector magnetograms. 
We conclude therefore, that non-linear force-free extrapolation meth- 
ods may not be used directly on observed vector magnetograms, in 
particular not on the very noisy transverse photospheric magnetic field 
measurements. The large noise in the transverse components of the 
photospheric field vector (~ the horizontal B x and B y at the bot- 
tom boundary) gives us the freedom to adjust these data within the 
noise level. We use this freedom to drive the data towards being more 
consistent with Aly's force-free and torque-free conditions. 

As a measure, how well a photospheric magnetic field agrees with 
Aly's criteria, the observed data and the smoothness condition, we 
define the following functional which adds up the chi-square deviations 
from all individual constraints: 



where 
Li = 

Lo = 



L 3 
L 4 



L = + H2L2 + ^3^3 + ^4^4 (6) 

(E B,B.) 2 + (E B v B ) 2 + (E Bl -Bl- B^j 2 
(j> (Bl -Bl- B 2 )J + {^y(Bl -Bl- B^ 
+ (^2yB x B z -xB y B^j 

E (B x - B xobs ) 2 + E (By - B yobs ) 2 + E (B z - B zobs y 
, p p p 

£ (^B x f + (AB y ) 2 + (AB 2 ) 2 j (7) 



The surface integrals are here replaced by a summation J2 P over a H 
grid nodes p of the bottom surface grid and the differentiation in 
the smoothing term is achieved by the usual 5-point stencil for the 
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2D-Laplace operator. Each constraint L n is weighted by a yet undeter- 
mined factor /j, n . The first term (n=l) corresponds to the force-balance 
condition, the next (n=2) to the torque-free condition. The following 
term (n=3) ensures that the optimized boundary condition agrees with 
the measured photospheric data and the last terms (n=4) controls the 
smoothing. The 2D-Laplace operator is designated by A. 

The aim of our preprocessing procedure is to minimize L so that all 
terms L n if possible are made small simultaneously. This will yield a 
surface magnetic field 



Besides a dependence on the observed magnetogram, the solution (8) 
now also depends on the coefficients fi n . These coefficients are for once 
a formal necessity because the terms L n physically represent different 
quantities with different units. By means of these coefficients, however, 
we can also give more or less weight to the individual terms in the case 
where a reduction in one term contradicts the reduction in another. 
This competition obviously exists between the observation term n=3 
and the smoothing term n=4. 

The smoothing is performed consistently for all three magnetic field 
components. For this purpose we need the derivative of L with respect 
to each of the three field components at every node q of the bottom 
boundary grid. We have, however, take account of the fact that B z is 
measured with much higher accuracy than B x and B y . This is achieved 
by assuming the vertical component invariable compared to the hori- 
zontal components in all terms where mixed products of the vertical 
and horizontal field components occur, e.g., in the Aly constraints. The 
relevant functional derivative of L is therefore 



B min = argmin(L) 



(8) 




+ 2^3 (B x - B xobs ) q + 2/x 4 (A (AB x )) q 



(9) 
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dL 



d(By) q 



= 2^i (j2 B v B ^j ( B *)g ~ 2 (E ( s . 
- 2 M2 2^x(i? 2 2 -^_ jB 2^ ( ^ )(; 

+ (Y,yB x B z - xByB^j {xB z ) q 
+ 2/x 3 (B y - B yobs ) + 2^ (A (AB y )) 



(10) 



c2L 



= 2^(B z -B zobs ) + 2/i 4 (A (AB Z )) 



(11) 



The optimization is performed iteratively by a simple Newton scheme 
which replaces 



x)q 



(B x )g ~ H 



dL 



(By) 



yjq 



(By)g ~ H 



(B z ) q ^ (B z ) q -n 



d{B x ) q 

dL 
d(B y ) q 

dL 
d{B z ) q 



(12) 
(13) 



at every step. The convergence of this scheme towards a solution of (8) 
is obvious: L has to decrease monotonically at every step as long as 
(9-11) has a nonzero component. These terms, however, vanish only if 
an extremum of L is reached. Since L is fourth order in B, this may not 
necessarily be a global minimum, in rare cases if the step size is handled 
carelessly it may even be a local maximum. In practical calculations this 
should, however, not be a problem and from our experience we rapidly 
obtain a minimum B mui of L, once the parameters /x n are specified. 

What is left is a suitable recipe to choose ji n . We have four param- 
eters fi n to choose and we restrict our freedom a little bit by giving 
the same weight fj,i = fi2D 2 = to the Aly momentum and torque 
constraints where D is the edge length of our numerical square box. 
Accordingly we will combine L\ + L2/D 2 = L\2- We are then left with 
three parameters only two of which are independent since only the ratio 
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of the parameters really counts. Multiplying all fi n with a common 
factor does not change the minimum solution B of the preprocessed 
photospheric data. We use this fact and specify [X\2 to 1/B^ ve where 
.B ave is the average magnetic field magnitude in our magnetogram. 1 

With the parameters ^3 and ^4 we can now control the relative 
influence of the observed data and the smoothing, respectively, to each 
other and to the Aly constraints. For the proper selection of these 
parameters we proceed in a similar way as it is customary for regular- 
ization techniques of inversion problems (e.g. Hansen (2000)) For this 
purpose we visualize the different solutions B mm (//3, ^4) for different \x 
parameters in a phase space of log(Li2), \og{Lj) and log(L4). Since the 
L-terms depend analytically in B, the solutions B m i n span a continuous 
2D surface in this 3D phase space. 

A survey of the [i parameters for one of the example computations 
(Low and Lou solution with noise model I) described in the next section 
shows that log(L3(B mm )) and log(L4(B mm )) are almost entirely deter- 
mined by the ratio of ^3 to ^4 while log(Li2(B mm )) depends on the 
absolute magnitude on /X3 and ^4. It is therefore convenient to display 
our results in a projection of the 3D phase space along the log(Li2)-axis 
as in Fig. 1. In this projection, the surface collapses nearly to a unique 
curve. 

An obvious limiting value is obtained for ^ — > 00 which yields B rmn 
= B bs w hh L12 ~ 10 _1 , L3 = and L4 ~ 10~ 3 . For smaller values 
of us we obtain a smoothed solution and, depending on 114 also satisfy 
the Aly criteria much better than the original observational data. The 
price we have to pay is that L3 attains finite values, which however is 
tolerable as long as L3 does not exceed a noise value L3(B no ; se ) (dotted 
line in Fig. 1). This value in our test case is known from the amount 
of noise added. For actual observations, it has to be estimated. The 
intersection of our solution surface with the noise level then defines 
the optimal ratio ^3 to 114 (marked by a rhombus in Fig. 1. We choose 
both numbers small to enforce a good compliance with the Aly criteria, 
i.e., give much weight to Lyi- This way, Lyi is reduced by 6 orders of 
magnitude, B is conveniently smoothed and yet it does not deviate from 
the observations by more than the instrument error. The smoothing will 
somewhat broaden prominent magnetic flux structures. A careful choice 
of the preprocessing parameters (as described above) ensures that the 
magnetic flux magnitudes and the corresponding magnetic field topol- 
ogy (which might become very complex for multiple magnetic sources, 
see (Schrijver and Title, 2002)) are not affected by the preprocessing. 

1 Please note that this is equivalent to a normalization of the magnetic field with 
B avo and the length scale with D. With this normalization we get L\ — L2 = L12 = 
1. 
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Figure 1. Survey of solutions to (8) in a phase space spanned by \og{Li) (abscissa), 
log(Z/3) (ordinate) and log(Li2) (color code). Every symbol corresponds to one so- 
lution B m i n . The optimal parameters are marked with a rhombus. Here we have 
ii2 = 7.6 ■ 1CP 7 , L3 = 6.3 • 10 -6 and La — 5.8 ■ 10 -5 . The optimal parameters are 
fj, 3 = 0.01 and ^4 = 0.005. 

4. Tests with the Low and Lou solution 

For a first test we use the semi-analytic solution found by Low and Lou 
(1990). This solution has been designed in particular to test non-linear 
force-free extrapolation codes. The solution contains free parameters 
and we use <3? = tt/2 and I = 0.3 (see Low and Lou (1990) for de- 
tails). We extract the bottom boundary of this equilibrium and use 
it as input for our extrapolation code (see Wiegelmann (2004)). This 
artificial vector magnetogram (see first row of Fig. 2) extrapolated 
from a semi-analytic solution is of course in perfect agreement with the 
assumption of a force-free field above (Aly-criteria) and the result of our 
extrapolation code showed a reasonable agreement with the original. 
Real measured vector magnetograms are not so ideal (and smooth) of 
course and we simulate this effect by adding noise to the Low and Lou 
magnetogram. We add noise to this ideal solution in the form 
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Figure 2. The first row shows the original vector magnetogram deduced from the 
Low and Lou solution. In the second row we added noise (model I) to the solution 
to simulate measurement errors. In the bottom row we applied our preprocessing 
routine to the noisy data taken from the center row. 
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Figure 3. a) Some field lines for the original Low and Lou solution, b) Potential 
field reconstruction, c) Non-linear force-free reconstruction from noisy data (noise 
model I) without preprocessing, d) Non-linear force-free reconstruction from noisy 
data (noise model I) after preprocessing the vector magnetogram with our newly 
developed program, e) Non-linear force-free reconstruction from noisy data (noise 
model II) without preprocessing, f) Non-linear force-free reconstruction from noisy 
data (noise model II) after preprocessing, g) Non-linear force-free reconstruction 
from noisy data (noise model III) without preprocessing, h) Non-linear force-free 
reconstruction from noisy data (noise model III) after preprocessing the vector 
magnetogram with our newly developed program. The color coding shows the normal 
magnetic field component on the photosphere in all panel. One can see, that after 
preprocessing (panel d,f,h) the agreement with the original (a) is much better than 
without preprocessing (c,e,g). 
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Noise model I: 

5E>i = nyr n - \/ r B~i, where n\ is the noise level and r n a random number 
in the range —1 ... 1. The noiselevel was choosen n\ = 5.0 for the 
transversal magnetic field (B x , B y ) and n\ = 0.25 for B z . This mimics 
a real magnetogram (see center row of Fig. 2) which has significantly- 
higher noise in the transversal components of the magnetic field. 
Noise model II: 

5Bi = ni-r n , where n\ is the noise level and r n a random number in the 
range —1 ... 1. The noiselevel was choosen ni = 10.0 for the transversal 
magnetic field (B x , B y ) and n\ = 0.5 for B z . This noise model adds an 
additive noise, independent from the magnetic field strength. 
Noise model III: 

SB 

SB Z = constant, 5B t = , „ t ° li ° — where we choose the constant noise 

V-^+^tmin 

level of B z to 1 and the minimum detection level SB tm i n = 20. This 
noise model mimics the effect that the transversal noise level is higher 
in regions with a low magnetic field strength, which is an immediate 
consequence of Eq. 4. 

The bottom row of Fig. 2 shows the preprocessed vector magne- 
togram (for noise model I) after applying our procedure. The aim of 
the preprocessing is to use the resulting magnetogram as input for 
a non-linear force-free magnetic field extrapolation. Fig. 3 shows in 
panel a) the original Low and Lou solution and in panel b) a corre- 
sponding potential field reconstruction. In Fig. 3 we present only the 
central region of the whole magnetogram (marked with black rectan- 
gular box in Fig. 2 because the surrounding magnetogram is used as 
a boundary layer (16 grid points) for our non-linear force- free code 
(Our non-linear force-free code is explained in detail in Wiegelmann 
(2004)). The computation was done on a 96 x 96 x 80 grid including 
a 16 pixel boundary layer towards the lateral and top boundary of the 
computational box.) In the remaining panels of Fig. 3 we demonstrate 
the effect of the different noise models on the reconstruction. The noise 
levels were chosen so that the mean noise was similar for all three noise 
models. Fig. 3 c) shows a non- linear force-free reconstruction with noisy 
data (noise model I, magnetogram shown in the center panel of Fig. 
2) and Fig. 3 d) presents a non-linear force-free reconstruction after 
preprocessing (Magnetogram shown in the bottom panel of Fig. 2). 
After preprocessing(Fig. 3 d) we get a much better agreement with 
the original solution (Fig. 3 a). Fig. 3 e) and f) show a non-linear 
force-free reconstruction for noise model II without (e) and after (f) 
preprocessing of the noisy data. Fig. 3 g) and h) show a non-linear 
force-free reconstruction for noise model III without (g) and after (h) 
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preprocessing of the noisy data. Similar as for noise model I we find 
that the preprocessed data agree better with the original Fig. 3 a). We 
check the correlation of the original solution with our reconstruction 
with help of the vector correlation function VC defined as 

VC = ^ V i ' Wi (14) 
^ 1 ^ ,/V.I Wi l 2 



Vi correspond to a reference field (Low and Lou solution) and Wi to 
the non-linear force-free field reconstructed with our code. 

Case remark preprocessing vector correlation 

a) reference 1.00 

b) potential 0.85 

c) Noise model I No 0.94 

d) Noise model I Yes 0.98 

e) Noise model II No 0.95 

f) Noise model II Yes 0.98 

g) Noise model III No 0.93 

h) Noise model III Yes 0.98 

The table confirms the visual inspection of Fig. 3. The correlation 
of the reconstructed magnetic field with the original becomes signif- 
icantly improved after preprocessing of the data for all noise models. 
We knew already from previous studies (Wiegelmann and Neukirch, 
2003; Wiegelmann, 2004) that noise and inconsistencies in vector mag- 
netograms have negative influence on the non-linear force-free recon- 
struction and the preprocessing routine described in this work tells us 
how to overcome these difficulties. 



5. Application to data from the Solar Flare Telescope (SFT) 

In this section we apply our method to vector magnetogram data 
of AR 7321 taken with SFT at the National Astronomical Observa- 
tory (NAO) in Tokyo. The SFT instrument is described in (Sakurai 
et al., 1995) and the evolution of the photospheric magnetic field in 
AR 7321 has been studied by Li et al. (2000). The top panel of Fig. 
4 shows the original data (original resolution reduced by a factor of 
two.) and the bottom panel after preprocessing. The structure of all 
components of the magnetic field is similar to the original. The mag- 
netogram is almost flux-balanced (efl ux = 0.045), but Aly's criteria 
are not fulfilled in the original (e^fy = 1-50). After preprocessing we 
get e ai v = 0.029 and the data are also somewhat smoother than the 
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Figure 4- Top panel: vector magnetogram of AR 7321 taken with SFT at October, 
26 1992. Bottom panel: After applying our preprocessing procedure. The resolution 
of the original SFT-magnetogram has been reduced by a factor of two to speed up 
the non-linear force-free magnetic field reconstruction. An inspection of the top and 
bottom panel shows that local small scale structures with strong flux density in the 
original lead to a broader structure with weaker flux density after preprocessing. 
Sharp boundaries between positive and negative flux regions are also broadened due 
to smoothing during preprocessing. 
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Figure 5. AR 7321 measured with SFT-data from October, 26 1992. a) Potential field 
reconstruction. b)Non-linear force-free reconstruction. The field lines start from the 
same footpoints within regions B z > in both panels. 



original. The L-values of the original data and after preprocessing are: 







original 


preprocessed 




(Aly-criteria) 


0.58 


3.6- 10~ 5 


Li 


(Data) 


0.0 


1.0- 10~ 5 


Lv 


(Smoothing) 


2.1 • 10~ 4 


5.2 • 10~ 6 



Fig. 5 a) shows a potential field reconstruction and Fig. 5 b) a 
non-linear force-free reconstruction based on the preprocessed vector 
magnetogram. The non-linear force-free computation was done in a 
96 x 96 x 80 box, including a boundary layer of 16 grid points towards 
the lateral and top boundary. Fig. 5 shows only the center cube of 
64 x 64 x 64 pixel which corresponds to the area marked with a black 
rectangular box in Fig. 4. 
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6. Conclusions 

Within this work we presented a method for the preprocessing of vector 
magnetogram data with the aim to use the result of the preprocessing as 
input for a non-linear force-free magnetic field extrapolation with help 
of an optimization code method. As a first test of the method we use the 
Low and Lou solution with overlaid noise with different noise models. 
A direct use of the noisy photospheric data for a non-linear force-free 
extrapolation showed no good agreement with the original Low and Lou 
solution, but after applying our newly developed preprocessing method 
we got a reasonable agreement with the original. The preprocessing 
method uses the high noise level in the transversal vector magnetogram 
components to drive the magnetogram towards boundary conditions 
which are consistent with the assumption of a force-free field above. To 
do so we use an minimization principle. On the one hand we control that 
the final boundary data are as close as possible (within the noise level) 
to the original measured data and on the other hand the data are forced 
to fulfill the Aly-criteria and be sufficiently smooth. Smoothness of the 
boundary data is required for the non-linear force-free extrapolation 
code, but also physically motivated because the magnetic field at the 
basis of the corona should be smoother than on the photosphere, where 
it is measured. 

Let us remark that non-linear force-free extrapolations are only valid 
as long as the plasma (3 is low and non-magnetic forces like pressure 
gradients and gravity can be neglected. This is usually the case in 
the corona (low j3 plasma, (3 w 10~ 4 ), but not in the photosphere 
(P ~ 1). A fully understanding on how the magnetic field evolves from 
the (non force-free) photosphere through the chromosphere and tran- 
sition region into the (force-free) corona would require more advanced 
models, at least magnetohydrostatic equilibria. Measurements of the 
photospheric magnetic field vector are not sufficient for such models 
and one needs additional observations regarding plasma quantities (e.g. 
density, temperature, pressure). A magnetohydrostatic approach is also 
required to model prominences where the magnetic field counteracts the 
gravitational force of the material they support. 

We apply our newly developed preprocessing program to data taken 
with the Solar Flare Telescope and use the preprocessed boundary data 
for a non-linear force-free field extrapolation. 
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